Dark energy and particle mixing 
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We show that the vacuum condensate due to particle mixing is responsible of a dynamically 
evolving dark energy. In particular, we show that values of the adiabatic index close to —1 for 
vacuum condensates of neutrinos and quarks imply, at the present epoch, contributions to the 
vacuum energy compatible with the estimated upper bound on the dark energy. 
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I. INTRODUCTION 
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^rfH ■ The experimental achievements proving neutrino oscillations [l|, 0] seem to indicate a promising path beyond the 
Standard Model of electro-weak interaction for elementary particles. On the other hand, an increasing bulk of data 
■ has been accumulated in the last few years paving the way to the emergence of a new standard cosmological model 
\ usually referred to as the concordance model. The Hubble diagram of Type la Supernovae (SNela), measured by both 
. the Supernova Cosmology Project [sj and the High-z Team [4] up to redshift z ~ 1, was the first evidence that the 
universe is undergoing a phase of accelerated expansion. Balloon born experiments, such as BOOMERanG ^] and 
, MAXIMA Q , determined the location of the first and second peak in the anisotropy spectrum of cosmic microwave 
i' background radiation (CMBR) pointing out that the geometry of the universe is spatially flat. If combined with 
O , constraints coming from galaxy clusters on the matter density parameter JIm, these data indicate that the universe 
, is dominated by a non-clustered fluid with negative pressure, generically referred to as dark energy, which is able to 
drive the accelerated expansion. This picture has been further strengthened by the more precise measurements of the 
^— I CMBR spectrum, due to the WMAP experiment , and by the extension of the SNela Hubble diagram to redshifts 
^ higher than 1 Several models trying to explain this phenomenon have been presented; the simplest explanation is 
, claiming for the well known cosmological constant A [9J. Although the best flt to most of the available astrophysical 
data Q, the ACDM model fails in explaining why the inferred value of A is so tiny (123 orders of magnitude lower) 
compared to the typical vacuum energy values predicted by particle physics and why its energy density is today 
comparable to the matter density (the so called coincidence problem). 

In this paper we study the possibility that a link between high energy physics and cosmology might be found in the 
mechanism of particle mixing. In our discussion we resort to previous investigations which led us to the conclusion 
that neutrino mixing might contribute to the dark energy budget of the universe [lol . [ll| . We show that the vacuum 
condensate due to particle mixing is responsible of a dynamically evolving dark energy. In particular, we show that 
values of the adiabatic index close to —1, both for vacuum condensates of neutrinos and quarks imply, at the present 
epoch, contributions to the vacuum energy compatible with the observed cosmological constant. We compute such 
jj] ' a value and show that the condensate could give rise also to the dark matter component of the Universe, besides 
pi , the accelerating one. Our discussion and conclusions rest on the OFT formalism for particle mixing, which has been 
extensively discussed in recent years in the literature [l^, [H, [ij, [H, [l^ [13, [iBl- For the reader convenience we 
summarize it in the Appendix A. 

The fact that the mixing phenomenon might be a source for the dark energy appears to be relevant from a genuine 
experimental point of view since, up to now, none of the exotic candidates for dark matter and dark energy, has been 
detected at a fundamental level. 

The layout of the paper is the following. In Section II we present the particle mixing condensate in the early and 
in the present epoch. In Section III we compute the fermion mixing contributions to the dark energy at the present 
epoch. Conclusions are drawn in Section IV. We outline the QFT formalism for fermion mixing in the Appendix A. 
In the Appendix B are reported useful computations. 



II. PARTICLE MIXING AND DARK ENERGY 



As mentioned above, experimental data indicate that the today observed universe can be described as an accelerating 
Hubble fluid where the contribution of dark energy component to the total matter-energy density is JIa — 0.7. 
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Moreover, the cosmic flow is "today" accelerating while it was not so at intermediate redshift z (e.g. 1 < z < 10) 
where large scale structures have supposed to be clustered. Thus, physically motivated cosmological models should 
undergo, at least, three phases: an early accelerated inflationary phase, an intermediate standard matter dominated 
(decelerated) phase and a final, today observed, dark energy dominated (accelerated) phase. This means that we 
have to take into account some form of dark energy which evolves from early epochs inducing the today observed 
acceleration. 

In this Section we show that the energy density due to the vacuum condensate arising from particle mixing can 
be interpreted as an evolving dark energy. The calculation here presented is performed for Dirac fermion fields in a 
Minkowski space-time. It can be extended to curved space-times, as it will be shown in a forthcoming work. 

Let us calculate the contributions p™" and p™" of the particle mixing to the vacuum energy density and to the 
vacuum pressure. Such a contributions are given respectively by the (0,0) and (j,j) components of the energy- 
momentum tensor of the condensed particles given in Eqs. (|MT|) - (|X23)) in Appendix A. 

The energy-momentum tensor density l'^i,{x) for the fermion fields ipi, i = 1, 2, 3 [20| . is 

: T^i^A^) — \ ■■ (^m(a;)7p 8^ *m(a;)) : (1) 

where = (V'lj "02 7 V's)"^ ^^id the normal ordering is with respect to the vacuum |0)m for the massive fields. Then 
the energy momentum tensor density of the vacuum condensate is given by 

^;r'(2;)=/(0(t)|:V(x):|0(i))/, (2) 
where 0(t)) / is the vacuum for the flavor fields (see Appendix A). 



A. Early universe epochs 
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In the early universe epochs, when the breaking of the Lorentz invariance of the vacuum is not negligible, p" 
presents also space-time dependent condensate contributions. This implies that the contribution p™" of the particle 
mixing to the vacuum energy density is given by computing the expectation value of the (0,0) component of the 
energy-momentum tensor : Too ■= J (Px : Tqo{x) : in the physical vacuum |0(t))j: 

^:^'7oo/(0(t)|:T°«(0):|0(0)/. (3) 

: Too ■ is for definition the Hamiltonian : H : in Eq. (jA24ll that, in terms of the annihilation and creation operators of 
tpi, ijj2 and ip3, is 

■■Too--=Y.T. [ ^'kc... + • (4) 

i r 

The notation is the one introduced in the Appendix (Eq. (jA3P ). Note that Tqo is time independent, moreover, within 
the QFT mixing formalism we have 

/(0|:Too:|0)/ = /(0(i)|:Too:|0(i))/ (5) 

for any t. We then obtain 

7^'^fe,4/(o|<.«i;,jo)/-(-/(o|/3;:|,/3L.|o); 

which, introducing the cut-off becomes 

- / dkk^ Wfe,l(s?2C?3ln2l' + S?3l^l3l')+Wfe,2 |-Sl2C23+e*^Ci2S23Sl3| 1^2 1' + 523^ l"^2! 

Jo '- ^ 

W/c,3 (|-Cl2S23 + e"* Sl2C23Sl3|^ |^23P + hl2S23 + Ci2C23Sl3|^ iWsP) ■ (6) 



mix 
rvac 
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Here ujk^i = ^JW^i^rril and the notation is the one introduced in the Appendix A for the CKM matrix elements (Eq. 
(jA2p ) and for the Bogoliubov coefficients (Eqs. (jAlSp and Eqs. (|A19|) 1. In any epoch, the energy density induced 
by the particle mixing condensate can be expressed as 

mix rptnix i ^rmix /y\ 

rvac vac ' vac \ ) 
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where the kinetic term T™'^ and the potential term V^l'^ are respectively given by 



vac ""^^^ ^^^^ ^^^^^^^^^^ ^^^^^^ 'vac 

2 r 1,2 1.2 



I" Jo 



Ai^ / 2 2 

+ I-C12S23 +e"^Sl2C23Sl3| |^23P+|si2S23+e'*Ci2C23Sl3| IW3P 



(8) 



and 



- - r&fc2[^(.?2C?3ln"2p+S?3l^l'3P) + — (|-M+e^'ci2W 
2 

+ ^ f I-C12S23 + e'"* Si2C23Sl3|^ IF23P + |S12S23 + 012023513^ l^iaP 
Wfc,3 V 



K2? + 4zclz\vh? 



(9) 



Eqs.® and ([1]) are obtained from Eq.® by using the relation ujk i = — h 

In a similar way, the contribution of particle mixing to the vacuum pressure is given by the expectation value 
of : Tjj : (where no summation on the index j is intended) in the vacuum \0{t)) f. 

P7:c^-lpVnfm\:T":\0{t))f. (10) 

Being 

■ - E E / '^'k ^ {Oh + , (11) 

in the case of the isotropy of the momenta we have T^^ — T^^ — T^'^, then 
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iCc = 17- / dkk' 



K 



— {slA, \V^,? + Sf, \Vt,n + — (I-.12C23 + e^' Ci2.23Sl3r \Vt,\' + 44 I ^2^3 P 
i^k,l ^k,2 ^ 



U2 
t^fc.3 



-C12S23 + e*'' Si2C23Sl3|^ |^23P + 1512523 + e"^Ci2C23Si3|^ |Vl3p 



(12) 



From Eqs.® and (HH), we define the adiabatic index ui™*^ = P™ac I Pm 'c ■ The plot of w™*^ as function of the 
momentum cut-off K (Fig.l) shows that = 1/3 when the cut-off is chosen to be ^ m where to is the largest 

of mi, TO2, TO3 and w™*^ goes to zero for K < ^toiTO2TO3 . 

This means that the condensate " mimics" the behavior of a perfect fluid of dust and radiation at the extreme values 
of the cut-off. From a dynamical point of view, it behaves as radiation in the relativistic regime (u;™" ~ 1/3) and 
as dark matter in the non- relativistic regime (ly™" ~ 0). Thus, in the early Universe and in the regions in which 
the breaking of Lorentz invariance of the vacuum is not negligible, the condensate could give rise to the dark matter 
component of the Universe. 

We note that according to this result, at the early universe epoch, the particle mixing condensate does not give 
contributions to the "standard" dark energy (the adiabatic index assumes, as we said, values in the range 

< w™" < 1/3). 

This gives the possibility to achieve the large scale structure formation as requested in a standard matter-radiation 
dominated regime and is in complete agreement with the WMAP results [2l| . Indeed, microwave light seen by WMAP 
from when the universe was only 380.000 years old, shows that, at that time, neutrinos made up 10% of the universe, 
atoms 12%, dark matter 63%, photons 15%, and dark energy was negligible. In contrast, estimates from WMAP 
data show the current universe consists of 4.6% of atoms, 23% dark matter, 72% dark energy and less than 1 percent 
neutrinos. 

The values of p™" and p™" which we obtain are time-independent since, as said, we are taking into account the 
Minkowski metric. Considering a curved space-time, time-dependence has to be taken into account but the essence 
of the result is expected to be the same (work is in progress on such an issue). 



B. Universe at present epoch 

At the present epoch, the breaking of the Lorentz invariance of the vacuum is very small and then p^'if comes 
almost completely from space-time independent condensate contributions (i.e. the contributions to the energy density 
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Figure 1: The adiabatic index lo™" as a function of cut-off K. 



of the vacuum |0)/ for mixed fields carrying a non- vanishing 9^ ~ fc^ = {oJk,kj) can be neglected). Then, in a flat 
space-time the kinetic term is negligible with respect to the potential ones : T™'^^ <C V^*^ and the 

energy-momentum density tensor of the vacuum condensate is approximatively given by 

r;r' ^ V,. E / ■■ M^)M^) ■ lo>/ = v,. pT- (i3) 

Since in a homogeneous and isotropic universe, the energy-momentum tensor density of the vacuum condensate can 
be written as T^^^^'^ — diag{p''J^'-'-^ jp"^^^ ,p^^^ ^p"]^^^ ), by equating this expression with Eq. (fT3|) and using 77^^ — 
diag{l, —1, —1, —1), we obtain the state equation: p™" ~ — Pa"j consistently with the vacuum Lorentz invariance. 

This means that the vacuum condensate, cor ning from particle mixing, contributes today to the dynamics of the 
universe with a cosmological constant behavior [llj. p™" computed from Eq. (fT5|) thus turns out to be 




which can be written as 




where eg = cos S. Note that p™" also depends on the CP violating phase i5. 

We observe that the value of the integral is conditioned by the appearance in the integrand of the \V^\'^ factors. 
The integral, and thus p™"' would be zero for \V^\'^ = for any |k|, as it is in the quantum mechanical (Pontecorvo) 
formalism [H, [13, [13, UHl ■ the present QFT formalism the jV^^p's account for the vacuum condensate (Eqs. (jA2ip - 
|A^3D ) and |V^^P goes to zero only for large momenta, getting its maximum value for |k| « ^/WiifrTj for any i,j = 1, 2, 3 

Proceeding in our calculation, we obtain that the integral diverges in K/^ as mf log {2Kf^/mj), with i, j = 1, 2, 3 
(see Appendix B). One also sees that '^'''^^x^'^^ c>c > for large Kj^. An interesting question to ask is how the 
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result pY^^ cx to^ \og{2Ki^/'mj), directly obtained in our approach, is related to the conjecture [265 that the small 
value of the cosmological constant oc {l{)~^eV)'^ is associated with the vacuum in a theory which has a fundamental 
mass scale m ~ 



III. PARTICLE MIXING CONDENSATE CONTRIBUTIONS AT THE PRESENT EPOCH 

In this Section we find a constraint on the cut-off on the momenta, at the present epoch, and we derive an expression 
of the adiabatic index of the particle condensates, w™" = p™"//°a"i ^ function of the cut-off. Then we show that 
values of the adiabatic index close to —1, both for vacuum condensates of neutrinos and quarks (denoted respectively 
with ix;^""'*'^ and w^"™") imply contributions to the vacuum energy p^"™" and p^"™" that are compatible with 
the estimated upper bound on the dark energy. 

The constraint on the cut-off is imposed by the very small breaking of the Lorentz invariance of the flavor vacuum 
at the present epoch. Indeed, by solving numerically the equations for T™" and V^'", given respectively by Eqs.® 
and we find that, in order to satisfy the condition T™-'^ <C t^", due to the very small breaking of the Lorentz 
invariance, the cut-off on the momenta at the present epoch must be chosen such that 

K\ <^ -^miTO2TO3 ■ (16) 

In particular, the exact value of the adiabatic index of the vacuum mixing condensates (of neutrinos and quarks) at 
the present epoch, tells us how much must be smaller than ^mim2mz- 

In the order to derive an expression of the state equation of the vacuum mixing condensates as function of K/y, let 
us consider the adiabatic expansion of a sphere of volume V . Let p denote the pressure at which the sphere expands. 
The total energy, E = pV , is not conserved since the pressure does work. Assuming that temperature and number of 
particles are constant, according to the first law of thermodynamics, the work done by p must be equal to the change 
in the total energy: dE = —pdV. That is pdV + V dp = ~pdV, that can be written as 

d[ip + p)V]^0, (17) 

from which 

const 

P + P^-^- (18) 
Then the equation of state into the sphere can be written as 

'^^- = T^K&— = 7r—r^ (19) 

r V f V ^ 

where C is a new constant. Note that w ^ —\ \i the volume is very large {V oo), that is, in the bulk of the 
Universe, i.e. far from the Universe "boundaries" . In collapsed regions {V 0) we have w — > . 

Eqs. p?|) - P^ hold for any fluid contained in an expanding volume V, when entropy, temperature, number of 
particles and electrochemical potential are assumed constant and p « const. Considering then the flavor vacuum 
condensate at the present epoch, taking into account the conditions: T^'^ <C V^"*^, and p™" ~ — p™" ~ V^"'^, from 
Eqs.® and nil) we have respectively pX*" = ^T" + ^a"" - ^a"" and p'l" = - pX*"- Thus the kinetic term 
is approximatively given by 

T™" ~ , (20) 

which means that, at the present epoch, the expansion of the universe leads to a smaller and smaller flavor vacuum 
condensate kinetic term. By using Eq. (|20p . the state equation for the flavor vacuum mixing condensate can be written 
as 

^rnix ^ P\ _ ^ (21) 

~Pa 

Eq. pTjl shows that, since at the present epoch 0, then w™" —1. Moreover, since T™" and p™" are 

function of the cut-off on the momenta Kj^, then Ea. ipT]) gives an expression of as function of Kj^: w™" = 

w^^^{K\). We now estimate the contributions given to the dark energy by the particle mixing condensates for 
different values of w™'^ close to —1, both for neutrino and for quark mixing condensates. 



6 



A. Neutrino mixing condensate contribution 

Let in Eq. (jA2[) represents the flavor neutrino fields: ^'J — {i^e, v^,Vr) and ^'m denotes the neutrino fields with 
definite masses, mi, m2, m^: = (i/i, 1/2, J^s)- The experimental values of squared mass differences and mixing 
angles are respectively: I^^ml^ = 7.9 x IQ-^eV"^, Am^g = 2.3 x IQ-^eF^, sf^ = 0.31, s^g = 0.44, sfg = 0.009 In 
the normal hierarchy case: \m^ \ 3> |toi^2|, we consider values of the neutrino masses such that the experimental values 
of squared mass difference are satisfied, as for example: mi — 4.6 x lO^^eV^, m2 = 1 x lO^^e^, ma = 5 x 
Then the condition Eg. pop for neutrinos reads 

KA^l.2xlO-^eV . (22) 

In Table 1, we report the contribution of the neutrino mixing to the dark energy and the corresponding 

state equation for different cut-offs satisfying the condition ([22ll . 











1.2 X 10 




1.1 X 10-^^ 


1.6 X 10-46 


-0.85 


4 X IQ- 


^eV 


1.2 X 10-4^ 


3.5 X 10-49 


-0.97 


3 X IQ- 


^eV 


0.3 X lO-''^ 


5.8 X 10-5" 


-0.98 


4 X 10~ 


*eV 


1.6 X 10-^2 


6.1 X 10-56 


-0.99 


4 X IQ- 




1.6 X 10-" 


6.2 X 10-63 


-0.99 



Table 1: Values of p^"'"'"' and u-^-™" for for different cut-offs. 



The result we find is that contributions to the dark energy compatible with its estimated upper bound: ™™ ~ 
10-4'^G'ey4 are obtained for values of the adiabatic index w^"™" of the neutrino mixing dark energy component: 

- 0.98 < wX"™" < -0.97 . (23) 

Eq. (|23p is in agreement with the constraint on the equation of state of the dark energy given by the combination of 
WMAP and Supernova Legacy Survey (SNLS) data: w — — 0.967j^n and with the constraint given by combining 
WMAP, large-scale structure and supernova data: w = —1.08 ± 0.12 |28|. 

A value of ly^-™" < —0.98 leads to negligible contributions of pj^"'"''^. The results we found are dependent on the 
neutrino mass values one uses. 



B. Quark mixing condensate contribution 



The quark mixing is expressed as: 




where V = 



is the CKM matrix 29]. In such a case, in Eq. (IA2 

\Vtd Vt, Vtb 

For the values of the quark masses given in Ref.|29|, the condition Eq. fTB]) for quarks reads 

Ka < 120Mey. 



= (d, s, b) and *T 



(24) 



(d',s',6'). 



(25) 



In Table 2, we report the contribution of the quark mixing to the dark energy and the corresponding state 

equation for different cut-offs satisfying the condition (pS)) . 



Ka 


pr""(Gey4) 


T«-™"(Gey4) 




UOMeV 


5.1 X 10-"^ 


3.5 X 10-7 


-0.3 


IQMeV 


2 X 10-1" 


1.4 X 10-" 


-0.93 


BOOKeV 


1.5 X 10-" 


1.8 X 10-21 


-0.99 


30i^ey 


1.5 X 10-22 


1.8 X 10-2S 


-0.99 


0.3eV 


1.5 X 10-47 


1.8 X 10-63 


-1 
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Table 2: Values of p\ ""'^ and w\ '"'^ for for different cut-offs. 

From Table 2, we find that the exact Lorentz invariance of the quark mixing condensate = — 1 (^j'?^™" 

is 16 orders less than V^^™"), at the present epoch, leads to a contribution to the dark energy that is compatible 
with its estimated upper bound: = 1.5 x \^~^'' GeV^ . We remark that very small deviations from the value 

^q-ra%x _ _y ^^^^ ^^^^ Contributions of that are beyond the accepted upper bound of the dark energy. 

The computation of p™'^ turns out to be sensible to small variations in the values of the particle masses and of 
Am^. Our results are therefore dependent on the mass values one uses. 

It is our future plan to compare the present approach with the one of Ref . {s^j based on string models of D-particle 
foam. 

In conclusion, we have shown that under reasonable boundary conditions the vacuum condensate from particle 
mixing can provide contributions to the dark energy compatible with the observed value of the cosmological constant. 
At the present stage, the novelty of the mechanism here proposed in the study of dark energy and the remarkable 
improvement in the computed order of magnitude without needs of postulating (till now unobserved) exotic fields, 
reveals that the QFT particle mixing scenario provides an interesting approach to the dark energy problem. 



IV. CONCLUSIONS AND DISCUSSION 



We have shown that the energy density due to the vacuum condensate arising from the particle mixing can be 
interpreted as an evolving dark energy that at present epoch has a behavior and a value compatible with the observed 
cosmological constant. This value is obtained by imposing values of the adiabatic index close to —1, both for vacuum 
condensates of neutrinos and quarks. Our discussion has been limited to the case of Minkowski space-time. In 
a forthcoming paper we will present the explicit computation in curved space-time. There we will show that the 
mixing treatment here presented in the flat space-time is a good approximation in the present epoch of that in FRW 
space-time. 

A very short summary of the observational status of art can aid to clarify the frame for our considerations and results. 
As mentioned in the Introduction, the data accumulated in recent years indicate that the universe is dominated by a 
non-clustered fluid with negative pressure (the dark energy) able to drive the accelerated expansion. As a tentative 
solution to the inadequacy of the mentioned ACDM model, many authors have replaced the cosmological constant 
with a scalar field rolling down its potential and giving rise to models referred to as quintessence 31| . Even if successful 
in fitting the data, the quintessence approach to dark energy is still plagued by the coincidence problem since the dark 
energy and matter densities evolve differently and reach comparable values for a very limited portion of the universe 
evolution coinciding at present era. In this case, the coincidence problem is replaced with a fine-tuning problem. 
Moreover, it is not clear where this scalar field originates from, thus leaving a great uncertainty on the choice of the 
scalar field potential. The subtle and elusive nature of dark energy has led to look for completely different scenarios 
able to give a quintessential behavior without the need of exotic components. In this connection, it has been observed 
that the acceleration of the universe only claims for a negative pressure dominant component, but does not tell 
anything about the nature and the number of cosmic fluids filling the universe [3^ . This consideration suggests that 
it could be possible to explain the accelerated expansion by introducing a single cosmic fiuid with an equation of state 
causing it to act like dark matter at high densities (giving rise to clustered structures) and dark energy at low densities 
(then giving rise to accelerated behavior of cosmic fiuid). An attractive feature of these models, usually referred to 
as Unified Dark Energy (UDE) or Unified Dark Matter (UDM) models, is that such an approach naturally solves, 
at least phenomenologically, the coincidence problem. Some interesting examples are the generalized Chaplygin gas 



33 1 , the tachyon field |3J| and the condensate cosmology [35| . A different class of UDE models has been proposed 
36| where a single fluid is considered whose energy density scales with the redshift in such a way that the radiation 



dominated era, the matter dominated era and the accelerating phase can be naturally achieved. Actually, there is 
still a different way to face the problem of cosmic acceleration. It is possible that the observed acceleration is not the 
manifestation of another ingredient in the cosmic pie, but rather the flrst signal of a breakdown of our understanding 
of the laws of gravitation |37l . Examples of models comprising only the standard matter are provided by the 
Cardassian expansion (3^, the DGP gravity higher order gravity actions [4l|, non - vanishing torsion fleld [i^ . 
higher-order curvature invariants included in the gravity Lagrangian [43|, etc.. 

This abundance of models is from one hand the signal of the fact that we have a limited number of cosmological 
tests to discriminate among rival theories, and from the other hand, that a urgent degeneracy problem has to be 
faced. The fact that the vacuum condensate originated by particle mixing provides contributions to the dark energy 
compatible with today expected value could contribute towards a solution of such a problem from both experimental 
and theoretical viewpoints. 
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Appendix A: PARTICLE MIXING IN QUANTUM FIELD THEORY 



The main features of the QFT formalism for the fermion mixing are here summarized (see |17| for a detailed review). 
The Lagrangian density describing three Dirac fields with a mixed mass term is: 



C{x) = ^^{x){i 



(Al) 



where ^'J = {•iPa,''Pbt^c) are the fields with definite flavors, and M = is the mixed mass matrix. Among the 
various possible parameterizations of the mixing matrix for three fields, we work with CKM matrix of the form: 



/ 



C12C13 



S12C13 



-S12C23 - Ci2S23Sl3e*'' C12C23 - Sl2S23Sl3e"' 



S23C13 



id ' 



(A2) 



V S12S23 - Ci2C23Sl3e*'' -C12S23 - Si2C23Sl3e"' C23C13 



with 



cos 9ij and stj = sin 9ij, being 9ij the mixing angle between tpi^ipj, S is the CP violating phase and 



— (V'lj ^2, '^^3) are the fields with definite masses mi ^ m2 ^ m^: 



k.r 



1,2,3, 



(A3) 



with al^^it) = al^^e-'"^^-^^ fi^^it) = (31]^ e^'^"-^', and u;k,^ = y/k^ + ml The operators a^^^ and (31, ^, z = 1,2,3 , r = 
1,2, annihilate the vacuum state |0)„i = |0)i (g) |0)2 ^ \0)3- Q^kilO)™ = Pki\^)m — 0. The anticommutation relations 



are 



: |zyf (a;),t/j^^(y)|^_^^ = ^^(x - y)(5a/3(5y, with = 1, ...4, and |ak,,,aq^} Sk^SrsStf, [p^^,, P'^j} = ^^^q^rsSrj, 
with i,j = 1, 2, 3. All other anticommutators are zero. The orthonormality and completeness relations are: u^^iU^^ i — 



s- ft 



„'-t 



-k,i 



= 0, and Er("k,^^k!^ + « 



„'-t 



-k,i^ 



I) — 1. Using Eq. (|A2p . we diagonalize the 



quadratic form of Eq. ljAip . which then reduces to the Lagrangian for three Dirac fields, with masses mi, m2 and m3: 

£(2;) = ^m{x) Md) *™(x) , (A4) 

where = diag(mi,m2,m^). 

The mixing transformation can be written as V"(a;) = Gg^{t)^f{x)Gg{t), where {a,i) = (A, 1), (B, 2), (C, 3), and 
the generator is now 



Ge{t) = G2i{t) G^^{t) Gi2(t), 



where 



Gl2{t) 


= exp 


012 


G23{t) 


= exp 




Gl3{t) 


= exp 


013 



(fx {^\{x)'4)2{x) - ll)\{x)'4}l{x 

(fx (^V|(a;)'03(a;) - il^l{x)il)2{x 
(fx {i^\{x)i,-i{x)e-'^ - v4(a;)V'i(a;)e"' 



(A5) 

(A6) 
(A7) 
(A8) 



At finite volume, Ge{t) is an unitary operator, G'g^it) = G-g{t) — G\(t\ preserving the canonical anticommutation 
relations; Gg^{t) maps the Hilbert spaces for ipi, 4>2 and 4'3 fields Hm to the Hilbert spaces for fiavored fields Hf. 
Gg (t) : "Hm 1^ Ti-f- In particular, for the vacuum |0)m we have, at finite volume V: 



\m)f^Gg\t) |o). 



(A9) 



|0)/ is the vacuum for 7i/, which we will refer to as the flavor vacuum. In the infinite volume limit the flavor vacuum 
\0{t)) f turns out to be unitarily inequivalent to the vacuum for the massive neutrinos |0)m [l3 |. This can be proved 
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for any number of generations by using rigorous mathematical methods [l5|. The non-perturbative nature of the 
flavored vacuum for the mixed neutrinos is thus revealed. 

Due to the linearity of Gg{t), we can express the flavor annihilators, relative to the fields ipaix) at each time, as 
(we use (a,*) = (A, 1), (i?, 2), (C, 3)): 

PiAt) ^ G-\t) (3l,{t) Ge{t) . (AlO) 
The flavor fields can be expanded in the same bases as v^: 



The flavor annihilation operators in the reference frame such that k = (0, 0, |k|) are: 

<Ait) = C12C13 al,{t) + s,2Ci, {\U^,\ alAt) + e^\V,\\ PLIA*)) + (iC^i^al + elV^ial /^^IsW) 



(A12) 



"L.B W - (ci2C23 - e'' S12S23S13) alAi) " (S12C23 + e^' C12S23S13) (\U^2\ "L.iW " ^1V^^\ P^At)] 

+ S23C13 {\U^s\ <3it) + e1^23l /^^^W) ' (A13) 

elicit) = C23C13 al,{t) - (C12S23 + S12C23S13) {{U^^l ^^1) - \V^%\ 

+ (S12S23 - e^' C12C23S13) {\U^s\ - el^al /^l.iW) , (AM) 

P-uAt) = C12C13 /3:k.l(i) + S12C13 {\U^2\ P-uA^) " ^ln'^21 al],{t)) + 6^' S13 (|C/i'3l P^^K^it) " ^1 ^^1^3 1 ^3 



(A15) 



(C12C23 - e-'^*' S12S23S13) /3:k,2W - (S12C23 + e"*' C12S23S13) (|C/^2l Ck,lW + 1^1^21 <i W 
+ S23C13 {\U^3\ PL^At) - |V^23l «L!3W) , (A16) 

Ck.C W = C23C13 /3:k,3 - (C12S23 + e"^* S12C23S13) (|C/2''3l /31k,2(i) + ^1 ^^23 I «M W) 

+ (S12S23 - e-'' C12C23S13) (|C/?3l /31k,iW + filial <i(0) • (A17) 

These operators satisfy canonical (anti)commutation relations at equal times, t/*^ and are Bogoliubov coefficients 
defined as: 

iV^l^f'^^MlJ^y (^!M±^ |k| |k| \ ^ ^^^g^ 



2a;fe^i / V 2tjfcj / V(a;A;j+mj) (wfe^^ + m^) 



|C/5P + I^,^P = 1, (A20) 

where j = 1, 2, 3 and j > i. The numbers of particles condensed in the vacuum are different for fermions of different 
masses: 

< = f{oimtnm)f = fm\N^cm))f = ^i^ci, \vt2? + \vt,\'' , (A2i) 
= f{m\K:\m)f - /(owiA^^^'now)/ = hM+e^' 012^23^131' 1^1^21' +^2341^^23!'. (A22) 

A/? = /(O(t)KnO(t))/ = /(0(t)|A^^^3nO(i))/ = |-Ci2523+e'^Si2C235l3r 1^231'+ 1^12*23 +e'^C^^^^^ 

(A23) 
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Since the vacuum |0)m for the massive fields is unitarily inequivalent to the vacuum |0(t))/ for the mixed (flavored) 
fields at time t, for any t, two different normal orderings must be defined, respectively with respect to |0)mj as usual 
denoted by : ... :, and with respect to \0{t))f, denoted by :: ... :: . The Hamiltonian normal ordered with respect to 
the vacua |0),„ and |0(i)) / are given respectively by 

H~^{0\H\0)„,^H+ 2 f d^k{u;k,i + iOk,2 + u^k.s) ^Y.Y1 I ^'^0;^..^!^"^ + /^C^-^y > (A24) 

i r 

■.-.H:: = H - f{0{t)\H\0{t))f =H + 2 J d^k{ujks + ^k,2 + ^k^s) - 4: ^ d'k [^^,1 (s^^c^g jy^k^p + p) 



W/c,2 (|-Sl2C23 + e"^Ci2S23Sl3| l^l2p+S23C?3 1^23 P 
^/c,3 ( I-C12S23 + e'* Sl2C23Sl3|^ |^23P + |si2S23 + Ci2C23Si3|^ | V^i 3 ^ 



(A25) 



The state \0{t)) f is a condensate of massive particle-antiparticle pairs. Note that the difference of energy between 
\0{t)) f and |0)m represents the energy of the condensed neutrinos given in Eas. (|A21[) - (|A23l) 



.(O(t)l : H : m)f = /(0(t)|ff|0(i))/ - ™(0|H|0)„, = 41^3^ [^^^^ (.a^^fa |V\^2p + s?3 \Vi%\') 



^/c,2 (^|-Sl2C23 +e^''ci2S23Sl3|^ \V^2\^ + sl^cf^ 1^23^ 
Wfe,3 I -C12S23 + e S12C23S13 II/23I + Sl2S23+e C12C23S13 ll/^gl 



Appendix B: BEHAVIOR OF pt^"" FOR Kk > mi.ma.mg 
By solving Eq. p^ . we obtain the explicit expression for which for Kp^ ^ mi,m2,TO3 reduces to: 



^ mU2shml{ms-m,) rV^j-m^ ^ 2slA, rnl{rri2 - rm) f^mj-mU 

Svr L J ml - to? V toi / Jrnl - to? V toi / 



2 ~ '"i 

„2 1^2 o™ ™ I o™2n i„„ / 2ifA\ / 2 2 , „2 \ _2 / ^^i^A 



Si3 ("^3 - 2to3771i + 2toi) log - (S12C13 + Sis) "^1 log 

' TO3 J ^ ' \ nil 

' '2Ka \ I ^ ml f 25^30^3 to^(to3 -TO2) ^j.^^^^ / - 
Y^toI - to| ^ W2 

2 I-S12C23 + e'* ci2S23Si3|^ m\(m2 - toi) _j / 
— tanli ' — — 



+ .?2C?3 (-^ - 2m2TO, + 2to?) log ^ } + ^|^f23£liqi^|_ZM Cretan 

V TO2 / J 27r L 7^2 _ ^2 



— m\ ^ "^2 

I Z(5 

;23 + e 

(If)} 



+ I-S12C23 + e*"^ C12S23S13I (m^ - 2m2TOi + 2TO2) log 

■2 |si2S23 + e"'ci2C23Si3|^ m|(m3-TOi) ( \Jm\ - m\ 



TO3 r. 
27r I 



^TO§ — TO^ ^ "^3 



2 I-C12S23 + e"' S12C23S13I to^(to3-TO2)^ ^ -\ ( \f' 
+ — ■ , tanh ' 



^ nil — m\ ^ ^3 

/2Xa 



+ Si2S23+e* C12C23S13 (2TO3 - 2TO3mi + TOj) log 



(A26) 



2 2/2 r, r. 2\ 1 f'2'Kf^\ /, t |2 2 2 \ 2 1 f'^K\\ 

+ S23C13 ('^3 - 2m3m2 + 2TO2) log I — — 1 + (^|-si2C23 + e C12S23S13I + S23C13J TO2 log I I 



,x |2 I ,1 |2\ 2 1 / '^Kk 

-ci2S23 + e S12C23S13 + si2S23 + e C12C23S13 ^3 log 

II I y V "^3 

2 /2K\\ 1 

-C12S23 + e*"^ S12C23S13 1 ('Tiz - 2TO3TO2 + 2TO3) log f — — j I . (Bl) 
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This quantity diverges in Kf^ as rrif log {2K\/mj), with i.j = 1, 2, 3. 
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